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Abstract 

Within the framework of the Landau-Lifshitz-Gilbert equation, using 
permalloy parameters, we study the statics and dynamics of flat cir- 
cular magnetic nano-structures with an in-plane magnetic vortex con- 
figuration, putting particular emphasis on the vorticity of the mag- 
netic state and on the (perpendicular) polarisation of the vortex center, 
which may be shifted with respect to the center of the circle. These 
binary degrees of freedom can in principle be used to manipulate two 
independent bits of information. 

Studying switching processes induced by in-plane and out-of plane 
field pulses we find that it is possible to switch the vorticity of the 
magnetic dot on a time scale of 40 ps in strong enough and short 
enough perpendicular external field pulses {Bf^^ « 0.5 T, duration 
« 40 ps). But for realistically small values of the Gilbert damping, 
only the vorticity can be switched this fast, and it turns out that it 
is better to dismiss the center of the circle totally, concentrating on 
flat 'nano-rings' with an inner radius Ri and an outer radius i?2. On 
these 'nano-rings' the vortex state is more stable, and with respect to 
the switching of the vorticity these structures have similar properties 
as circular dots. 
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1 Introduction 



Recently there is a strong interest in using small magnetic structures as 
storage elements and for fast magnetoelectronic computation (MRAMs, 
imi). R. Cowburn, p|, has studied experimentally and by simulation 
arrays of flat magnetic structures with circular or square individual 
geometries, and M. Schneider et ai, 0, have considered more general 
geometries of the dots. The present state of the art concerning the 
static behaviour of magnetic nano-dots and their arrays is contained in 
the very recent papers of Ross et al, 0], Metlov and Guslienko, 0, and 
Guslienko et al, 0, 0. Concerning the dynamics we explicitly mention 
the even more recent papers of Gerrits et al., and Novosad et ai, 
10, and the references therein. 

If for a given thickness the radius R of the fiat dot is smaller than 
a critical value Rc (see below), then the dot is in an in-plane single- 
domain state, whereas for R > Rc, the dot is magnetically in a vortex 
state, see 0. 

Below, for fiat dots, i.e. for small values of their aspect ratio a := 
(th/R) -C 1, where th denotes the thickness of the dot, at first (for given 
R) the critical thickness tc is calculated for structures with permalloy 
material parameters, and a considerable stability range for the mag- 
netic vortex state is obtained even in rather tiny nano-structures. By 
application of a sufficiently small in-plane static magnetic field the 
magnetic state of the dot does not change topologically, i.e. the vortex 
center shifts, however the vorticity is still well-defined, and also the 
central direction of the magnetization (so-called " central polarisation" 
±) is unchanged. At the vortex center and in the region nearby, the 
magnetization (which must be constant in magnitude) is out-of-plane, 
whereas far from the center it is of course in-plane. 

Therefore the magnetic vortex state of the dot has a fourfold topo- 
logical degeneracy, since in principle one can measure (and flip) inde- 
pendently both the vorticity (r/l) of the state and the central polarisa- 
tion (+/— ) of the vortex; i.e. one has two bits per dot for a 'vortex-state 
dot', instead of only one bit for 'single- domain dots'. 

At present, the vorticity of the dot can be measured e.g. by electron 
microscopy in the Lorentz mode while the central polarisation can be 
measured e.g. by magnetic force microscopy, |T0[: Of course, for the 



applications one should find more convenient methods; but this task is 
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not considered here. In any case, the possibihty of storing and switching 
two independent bits of information instead of only one, and the fact 
that one does not rely on the extremely small single-domain dots, gives 
a strong additional motivation to study the statics and dynamics of 
circular dots in the vortex state very thoroughly. 

In the following chapter 2 we introduce the basic micromagnetic 
tools for such a study, the Landau-Lifshitz free energy functional and 



the dynamical Landau-Lifshitz-Gilbert equation, [|TT|. In chapter 3, 
the behaviour of flat circular dots under static external magnetic in- 
plane fields is considered, with particular emphasis on the transitions 
between the vortex state and the in-plane single-domain state, whereas 
in the next chapter 4 we consider the dynamics, both for in-plane fields 
and perpendicular field-pulses. The switching-dynamics is exemplified 
by some kind of "movie representation" for typical cases, which show 
that the time scale, r ~ 40 ps, for the precessional switching of the 
vorticity by strong perpendicular field pulses, H'^^^ ~ 5000 Oe over this 
duration, is extremely short, whereas the polarisation of the vortex 
core cannot be reliably switched at all on these time scales, see below. 
Therefore at the end, it is natural to skip the vortex center totally by 
considering in chapter 5 fiat nano-rings, i.e. circular structures with an 
inner radius Ri and an outer radius R2, which are both small, i.e. of 
the order of 50 to 150 nm, but large compared to the thickness th- We 
find that in this case the stability of the information is significantly 
improved and the switching dynamics is as fast as before. 



2 Basic equations 

Let us write the magnetic polarisation J(r, t) as J{r,t) = Jg ■ d{f,t), 
where Js = \J\ is the constant magnitude and d{r,t) the position- 
and time-dependent direction of the magnetic polarisation. We neglect 
in the following the small (in-plane) uniaxial magnetic anisotropy of 
Py, taking only the exchange interaction, the magnetic dipole-dipole 
interaction, and the Zeeman interaction (i.e. the interaction with the 
external magnetic field) into account. As a consequence, the Free En- 
ergy of the magnetic system K considered is obtained by minimization 
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of the following " Landau-Lifsliitz functional", see e.g. [|T2], |T3]: 



Fll:= J d\/-|A^(|^)2-J,-(l#°'P-a+i?-*-a)| 



(1) 



Here V{K) is the integration volume of the magnetic system K, A 
is the "exchange constant", Jg is the above-mentioned saturation polar- 
isation (we use A = 1.3 x 10~^ erg/cm and Js = 1 Tesla for Py, unless 
otherwise stated); H'^^^ is the external magnetic field and if^'P the in- 
ternal dipolar field (also called: "stray-field") which is generated by 
the magnetic moments of the system according to the magnetostatic 
equation H^^^{r,t) = —grad(j)m{r,t), with V^0m = ^ div ( • a) 
and the corresponding boundary and interface constraints. (Note that 
at this place, even on a ps time scale, we can safely use the magneto- 
static equations without taking into account retardation effects, since 
e.g. in the time of r = 1 ps a vacuum electromagnetic wave would 
traverse a distance of 300 /xm, which is far beyond the length scales of 
the structures considered, P^.) 

The dynamic equation of motion corresponding to (|T]) is 

^ = -j,dxH^\ (2) 

where 70 = 9^-^ is the " gyromagnetic ratio" (/i^ is Bohr's "magneton" , 
h Planck's constant divided by 2tt, and g the Lande factor, g = 2 for 
spin magnetism, i.e. in this case 70 has the value 1.176 ■ 10'' s~^/Oe), 
whereas H^^ is the space- and time-dependent effective magnetic field, 
around which the magnetic polarisation precesses, namely 

^eff ^ ^ext ^ ^Dip + ^ - ^ ■ ^ . (3) 

Js 7o 

Here - as already mentioned - ac is Gilbert's damping constant, 
|rT|; i.e. the last term in equation (|^) describes an isotropic velocity 



dependent phenomenological damping leading to a gradual decay of 
the precession amplitude. In principle, and particularly near interfaces, 
this damping should be anisotropic, see e.g. [T^, |TB[, but in the following 



this damping anisotropy is neglected, and also the thermal fluctuation 



4 



fields, which would be proportional to (T ■ ac) 2 , are neglected as usual 



(but see the recent papers of Brown et al. and Nowak et al, [|l7i P^). 
Also deviations from cubic symmetry of the magnetic material, which 
would lead to a more general density of the exchange energy, namely 

i Aj§--^,^Te neglected. 



3 Static behaviour of a magnetic dot: Vortex 
state versus in-plane single-domain state 

3.1 Vanishing external field 

At first we consider the static behaviour of a flat circular magnetic 
dot for vanishing external field; in particular we consider the phase 
boundary between the vortex state and an in-plane single-domain state, 
1l9| , minimizing the "Landau-Lifshitz Functional" (|I]) with the ansatz 



Oixi^) = — sin (y9(r) • (1 — a^(r)) 2 

ay{r) = +cos<^(r)- (l-a^(r))5 (4) 

Our special ansatz for the 'profile function' ^^(r) is an improved ver- 
sion of an approximation by Feldtkeller and Thomas, pO[)' 
following considerations: 

For r Im, where Im '■= {^Ai^)^ is the 'magnetic exchange length', 

s 

one expects the behaviour a;^(r) = 1 — const, -r^ — ... (where the dots de- 
note higher-order terms); beyond that range one will observe an expo- 
nential spatial dependence with decay length Im, while asymptotically 
for r ^ Im some kind of 'Mexican Hat' behaviour will apply, i.e. after 
the pronounced exponential decay azir) should become negative and 
approach the r-axis from negative values. However, this negative part 
should be proportional to the magnetic moment of the central part of 
the dot, i.e. ~ J^vr/^t/j, and therefore it vanishes cx for very fiat dots. 

(Our ansatz for az(r), see below, is more accurate than that of 
eq. (5) of Usov and Peschany, |]^, who neglect the magnetostatic en- 
ergy generated by the vortex core completely to obtain their simplified 
profile function, which contains a singularity of the derivative at the 
so-called core radius and is constant for larger r, i.e. it does not con- 
tain the tail nor the negative part of the profile. Our calculation is also 
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more accurate, see the results, than that of Gushenko and Metlov, 



and Gushenko et al, who take the magnetostatic field produced 
by the vortex core into account, but use basically the simplified core 
profile of ^ instead of the more correct one of below.) 

In any case, as a systematic approximation, expanding the ansatz 
of Feldtkeller and Thomas, [^], we use in eq. (§): 



^2 „2 



az{r) = c-e + (1 — c) ■ e (5) 

where c depends on the thickness of the sample and is obtained by an 
energy-fitting procedure leading to 

c ^ -0.188 + 0.708 ■ (1 + 0.2 • + 0.0278 ■ T^)'^ . (6) 

Here Th := is the reduced thickness of our structures. Note that 

I'm 

c ~ 0.5 for very fiat dots {th <C R). (Eqs. (|) and (H) are of course only 
approximately valid, since 'Mexican hat profiles', which one would get 
for thicker dots, are explicitly excluded for —0.188 < c < 0.52.) 

In fact, the gross behaviour of the profile function for our fiat struc- 
ture is given by the above-mentioned ansatz (|^): This is exemplified in 
Fig. 1, where the results of a numerical simulation of (Xz{r), performed 



with the OOMMF program, [^], for fiat circular permalloy structures 
with R = 175.5 nm and a thickness of th = Q nm, are compared with 
(1^) and d^). Details are given in Fig. la and in the corresponding fig- 
ure caption. In Fig. lb an effective radius i?efr of the central region of 
the vortex is presented as a function of the reduced thickness Th. This 
effective radius -Reff is defined by the equation R^g = f ■ — s-, taking 

into account the fact that the polar angle 6 of the magnetic polarisation 
is zero rsp. | in the vortex center rsp. far from it. With (§) one gets 
Rcs = ^^^g ^ ■ Im, which is plotted in Fig. lb. This means that i?eff is 
approximately between 1.5 and 3lm depending on the thickness th- 
Our ansatz az{r), eq. (^, for the above-mentioned profile function 
leads to analytical formulae for the exchange energy density f^^'^^ = 

3 

y I J2 (^)^dV, and we obtain for thin dots with c = 0.5 approxi- 

mately: /^^^^ ^ ^ ■ [2.67 + In ^g^]. 

For circular single-domain dots with homogeneous magnetization 
there is only the stray-field energy. For in-plane states, magnetic poles 
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arise exclusively at the side-face S (i.e. r = R), and we obtain for the 
stray-field energy density 

/Dip .= 1 . [ [ i^jr)- n{r)) ■ {a{r)-n{r)) ^2g^2g' 

V Sttuq J J \f—r"\ 

^ S S' ' ' 

By numerical evaluation, followed by a fit, we get: ^ ■ [2.41 ■ 
ln(l + 5.31(||)0-88)]. 

From these results one obtains for the critical thickness tc{R), where 
the energies of the vortex state and the in-plane single-domain states 
agree, with the ansatz (||): 

(8) 

Here the term ~ 0.5 comes from the dipolar energy of the central part 
of the vortex. 

In Fig. 2, this result for the critical thickness is plotted against 
R (the solid line) and compared with the results of an unpublished 
simplified calculation, (see |^ the dashed line, i.e. tcfi{R)), where 
the fact that in the center of the vortex is different from zero has 
been totally neglected (yielding a wrong asymptotics for R <C Im), 
and where the magnetostatic energy has been calculated by replacing 
the circular dot of radius R and thickness th by a prolate ellipsoid. 
It turns out that this simple calculation is not bad and captures the 
essentials, although there is still a significant discrepancy w.r.t. the 
direct numerical simulation e.g. by the OOMMF program. In contrast 
the better results with our ansatz (^ for the profile functions can hardly 
be distinguished, to the accuracy of the plot, from the results of the 
OOMMF simulation. 

(Additionally, Metlov and Guslienko, in their recent paper Q on 
thicker cylindrical dots, have determined the transition line between in- 
plane and out-of-plane single domain states, see Fig. 2 in [Q. However 
for our fiat dots with t/j <^ i? this addition is not important. In fact, 
in our Fig. 2, the additional phase boundary would be represented 
by a very steep line almost parallel to the abscissa axis, which would 
terminate at our 'solid line' phase boundary at a thickness of ~ 43 nm. 
This is far beyond the thicknesses we consider.) 
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3.2 In-plane external fields 

For in-plane external fields, the vortex deforms, and its center shifts 
by an amount A (in x-direction, if the external field is applied in the 
^/-direction; actually, see Fig. 8, below, the dynamic behaviour is more 
complicated). For simplicity, in the following approximation we only 
take into account the shift. Here both A < R and also A > R are con- 
sidered, see below). With this approximation we minimize the Landau- 
Lifshitz functional. For the 'stray-field energy density' we again have 
to evaluate (|^ with the (fictitious) magnetic poles at the side-face S; 
additionally there is now the Zeeman energy density 

/Zeeman 7|7jext| ^(i?. A, /„, c) 

/ =-Js\H I- , (9) 

where T{R, A, Im, c) is a function of the parameters R, A, and c. 
Note that - as already mentioned - A can be larger than R, although 
one expects that the transition from the vortex state to the in-plane 
single-domain state just happens in the region A ^ R. 

In Fig. 3, again for a circular Py dot of R = 150 nm and th = 4: nm in 
an external in-plane field of Hy^^ = 100 Oe, the energy densities (dipole 
energy, Zeeman energy, exchange energy, and total energy) are plotted 
as a function of the shift parameter A. The equilibrium position of the 
vortex center would correspond to the minimum of the total energy, 
which appears at A ^ 0.4 i?. With increasing A {—>■ R), in Fig. 3 
there is an energy barrier, separating the in-plane single-domain state 
(i.e. A ^ R) from the deformed vortex state (i.e. A < R). However 
this energy barrier depends on H^^^, and with increasing Hy^^ it decays 
to zero as shown in Fig. 4. 

Fig. 4 yields quantitative results for the 'annihilation fields' of the 
vortex state rsp. for the energy barrier separating the vortex state from 
the in-plane single-domain state, if the in-plane external field is not 
strong enough for the transition. Since the deformation of the vortex 
during the shift is neglected in our approximation, this energy bar- 
rier is overestimated. In any case, near 'rough boundaries' the actual 
behaviour will be too complex for a simple calculation, see ||25|| . 

In Fig. 5, we compare our calculations with experimental results 
for the annihilation fields of flat circular Py samples, plotted over the 
double reciprocal aspect ratio 2R/th, see [^. In our calculation we 
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have used the experimental values A = 1.3x 10~^ erg/ccm, Js = 0.9676 
Tesla, zero anisotropy, and two thicknesses (th = 15 nm and th = 8.3 
nm). As can be seen from the figure, in agreement with the above- 
mentioned arguments, our theoretical estimates are systematically ~ 
20 — 30 % too high, but a large part of the deviation is probably due 
to the fact that in the experiments a non-neglegible perpendicular field 



component was present, |2^; moreover, since the experiments are at 



room temperature, the infiuence of thermal fields may also have been 



non-neglegible, see ref. |16, E8[. In any case, in the paper of Guslienko 



and Metlov, [^, who also perform a comparison with the experiments 



of p6|, the theoretical results are significantly overestimated at small 
values of R/th, contrasted with our Fig. 5. The discrepancy is due 
to the fact that the authors of use basically the simplified profile 
function of Usov and Peschany, 0, which should not be applied for 
R < th- But globally, i.e. averaged over all values of R/th, in 0, the 
agreement between the calculated vortex annihilation fields and the 
experiments looks better than that one presented in our Fig. 5, in spite 
of the drastic overestimation at small values. 



4 Dynamics 

This topic is particularly emphasized in our paper; to use both the 
vorticity and the central polarisation of fiat circular dots for information 
processing one should be able to switch both properties reliably and 
independently, and - of course - fast, i.e. significantly below ns time 
scales. 

Unfortunately, the vorticity of the static magnetic state cannot be 
directly infiuenced experimentally with circular dots. However, M. 
Schneider et al. have shown how one can easily infiuence the vorticity 
of the magnetic state by using dots, which are not perfectly circular, 
cutting-off a segment of the structure, see This refers to static or 
quasi-static conditions. But in the following we consider fast dynamic 
changes, at first by time- dependent in-plane fields, and then by time- 
dependent perpendicular fields : 
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4.1 Time-dependent in-plane fields 

As usual, we consider a circular Py dot {A = 1.3 x 10~^ erg/ccm, Jg = 
1.08 T, neglegible anisotropy) with R = 150 nm and 4 nm thickness. 
Without external field, the dot is in a stable vortex state. At time t = 0, 
an in-plane external field Hy^^ is applied, and the 'critical duration' tcrit 
is determined, i.e. the time until the vortex has been shifted to the edge 
of the dot (A = i?), so that the vortex state is either 'lost' to an in-plane 
single-domain state, A > i?, or restored (A — > 0) after the field H^^^ 
is switched off. The results (which have been obtained by computer 
simulations with the OOMMF code, [^) are presented in Fig. 6; they 
depend on Gilbert's damping parameter ac in a rather subtle way 
(see below). At first, we see from this figure that to reach critical 
times below t = 1 ns one needs quite strong external fields, namely 
if^^* > 300 Oe to get below 800 ps, and > 500 Oe to get below 600 ps. 
Moreover, the information contained in the polarisation of the vortex 
center is only preserved for unrealistically large ac (e.g. for aa = 1 the 
vortex center is only stable for fields Hy^^ < 1100 Oe, i.e. up to the end- 
point of the lowest curve in Fig. 6.) For ac = 0.5 the corresponding 
(second-lowest) curve in Fig. 6 ends already at Hy^^ = 650 Oe with 
^crit = 500 ps, and for ac = 0.01 we have only obtained one point in 
Fig. 6, namely with tcrit = 1000 ps and H^^^ = 200 Oe. If the vortex 
state is 'lost' to an in-plane single-domain state by the application of 
the in-plane field Hy^^, one would need roughly at least another 'critical 
duration' to generate a new vortex state beginning from the boundary 
of the circular dot and propagating the core of the new vortex to the 
center; maybe this new vortex has inverted vorticity, inverted according 
to the systematic but slow prescription of p, |2^. In any case, this 
vorticity-switching would be too slow, and too complex. 

In fact, one sees from Fig. 7, where in an {ac, i^^^*)-plane the 
boundary line separating stable and unstable regions of the (deformed) 
vortex state is plotted, that for realistically small values of ac the 
vortex center is unstable for fields above 200 Oe, i.e. the independent 
information contained in the central polarisation of the vortex gets lost 
by the application of Hy^^: Stability of this 'polarisation information' 
can only be obtained for smaller in-plane fields Hy^^, which also means 
'critical times' larger than 1 ns. 

In the next figure. Fig. 8, which corresponds to some kind of ^ movie 
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sequence', the particular case of ac = 0.02 is considered in detail, for 
— 300 Oe. (The observation that the vortex core moves in the 
'north-east' direction, although from the fact that the switching-field is 
applied in the y-direction one would expect that it should move exactly 
'eastwards', comes from the following: In the first ~ 100 picoseconds 
the vortex core, in fact, starts moving to the right, but at the same 
time some kind of 180° domain wall is formed on the segment between 
X = and x — R; under the infiuence of the Hy-fie\d this domain wall 
tends to move upwards, and after 600 ps it can still be identified in Fig. 
8 to the r.h.s. of the vortex core. In any case, since the 'domain wall 
segment' is fixed to the (slower) vortex core, the "north-east" motion of 
the core results, as a consequence of the complicated coupling between 
all spins in course of the dynamics. Further details are only seen if the 
'movie' corresponding to Fig. 8 is observed on a shorter time-scale.) 

The total simulation time of Fig. 8 is 1.5 ns: One concludes that 
vorticity annihilation is interesting and complicated, but rather slow 
(t„it ~ 1.5 ns); particularly, the information contained in the central 
polarisation is lost by applying the field Hy^^ for the critical duration, 
and one is far from simply 'switching' the vorticity. (These conclusions 
would also apply for larger values of ao, although for aa — 1, which is 
unrealistically large, one would gain a factor of 3 in speed.) 

As a consequence, for vorticity-switching one should not use in- 
plane fields, but instead one should work with perpendicular fields, as 
discussed in the following. 

4.2 Time-dependent perpendiculeir fields 

We have finally performed simulations, starting from the planar equilib- 
rium vortex configuration, with pulsed external 'switching fields' which 
are now perpendicular to the plane. We consider the same structure as 
before. The applied field pulse was very strong, Bf^^ = HoHf^^ — 0.5 
T, i.e. i/f*=5000 Oe, but also very short (see below). The Gilbert 
damping was assumed to be as small as in Fig. 8, i.e. aa — 0.02, but 
the pulse time of Hf^* was only 42 ps, i.e. only a small fraction of the 
total simulation time of the corresponding 'movie sequence', which we 
calculated for that example. This time the vorticity is really inverted, 
and the switching is very fast, 42 ps, and without "ringing" after the 
end of the field-pulse. 
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In Fig. 9 we consider the effective 'switching time' tsw = ti, i.e. 
the time needed to turn around the spins - in a spatial average - 
by 180 degrees, systematically as a function of Hf^^. For vanishing 
damping this corresponds to a gyromagnetic precession by an angle 
71, i.e. ti = i^^i'^e^t • This is shown in Fig. 10, which presents results 
from the first part of our 'movie'. After a duration of ti , 42 ps in our 
case, the field is switched off; the system then 'wiggles' around the new 
equilibrium vortex state and gets slowly into equilibrium after several 
ns. But the 'wiggles' are small enough that the new equilibrium state 
can be clearly recognized already after 42 ps, although the system is still 
far from equilibrium, | |30| . 

(Note that in Fig. 10 the exchange energy does not change with 
time, in contrast to the magnetostatic energy, although after ~ 21 ps 
the transient spin configuration has changed from a vortex state to a 
'hedgehog state'; this invariance of the exchange energy is in agreement 



with a 'topological statement' of W. Doring, |31[| .) 

In Fig. 11, the 'wiggling' of the spatial average of az is presented. 
Finally in Fig. 12 we plot, what happens to the polarisation of the vor- 
tex center under the influence of a 'tilted vertical field pulse', involving 
a strong perpendicular component Hf^^ = —1000 Oe, applied between 
t = and 60 ps, accompanied by a simultaneous weaker in-plane pulse 
of Hy^^ = 150 Oe between t = and t = 30 ps. We find that at first the 
central polarisation is inverted from (+) to (-), but then it returns to 
the starting state (+). In fact, we find it difficult to control the central 
polarisation independently and with similar speed as the vorticity. 

That very strong vertical fields in the kOe range are necessary for 
polarity-switching of the vortex core is also found in a recent letter of 
Kikuchi et al, [g^. 

We have seen that the fast switching of the vorticity by vertical fields 
works successful, but the 'switching field pulses ' are in fact very large 
in strength, ^ 5000 Oe, and very short in duration, ^ 40 ps. However, 
concerning the independent switching of the polarisation of the center 
of the vortex, we always have found (for realistically small values of 
ac ) that the central polarisation - in contrast to the vorticity - can 
hardly be controlled even by strong field pulses oriented perpendicular 
to the plane; therefore in this case we are unsuccessful and do not plot 
the results of our simulations. 
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4.3 Vortex structure in flat circular dots: discussion 



Here we collect our results for the flat circular Permalloy (Py) dots 
considered in this section 4: 

1. The switching fields are quite large, i.e. necessary values are higher 
than 300 Oersteds (for in-plane fields), and 5000 Oersteds for the 
case of perpendicular fields. 

2. For realistic values of the Gilbert damping parameter ac see 



11], only the vorticity can be reliably switched, p3l, but not the 



central polarisation of the vortex core. 

The time scale necessary for the precessional switching of the 
vorticity is very short, r ~ 40 ps, for the case of perpendicular 
field pulses. In that case, after the switching the magnetization 
"wiggles" slightly, i.e. there are relatively small oscillations, below 
±10%, around the new equilibrium. These small precessional os- 
cillations decay rather slowly on a much larger ns time scale. But 
principally the information contained in the new state can already 
be obtained after 40 ps with sufficiently strong perpendicular field 
pulses without necessity to suppress the 'wiggling', . 

For the case of in-plane switching fields the field strengths for the 
dynamic annihilation of the vortex state are an order of magnitude 
lower, By^^ ~ 300 Oersteds, and they need not to be pulsed; but 
the time scales for vortex-switching with in-plane fields are also 



an order of magnitude slower, namely > 500 ps, |35 



5 Vortex structure in flat circular 'nano-rings' 

As a consequence, at the end, we dismiss the vortex center totally, al- 
though this implies the definite loss of the the possibility, to manipulate 
a second bit. We consider fiat circular 'nano-rings' with an inner radius 
Ri and an outer radius R2. In this case, see Fig. 13, with increasing _Ri 
the energy density of the vortex state decreases, while that of the ho- 
mogeneous in-plane state ('single- domain state') increases significantly, 
i.e. the stability of the vortex state is enhanced. 
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At the same time the switching of the vorticity by perpendicular 
field pulses is hardly influenced by the size of Ri. In Fig. 14 we present 
simulations for Py dots with R2 — 150 nm and th — 4: nm, comparing 

our former results for i?i = with similar 'ring results' for Ri = 50 
nm, as a function of the switching field Hf^^. Obviously the switching 
time ti is practically not changed. 

6 Conclusions 

Within the framework of the Landau-Lifshitz-Gilbert equation, using 
permalloy parameters, we have studied the statics and dynamics of 
fiat circular magnetic nano-structures with an in-plane magnetic vor- 
tex configuration putting particular emphasis on the (perpendicular) 
polarisation of the vortex center (which may be shifted with respect to 
the center of the circle), and on the (in-plane) vorticity of the magnetic 
state. Studying switching processes induced by in-plane and out-of- 
plane field pulses, we find that it is possible to switch the vorticity 
of the magnetic dot on an ultra-short time scale of 40 ps with strong 
enough and short enough perpendicular external field pulses (strength 
/loHf^^ pa 0.5 T; duration pa 40 ps); but for realistic values of the 
Gilbert damping constant ao, only the vorticity can be switched this 
fast, but not the magnetic polarisation = ±1 of the vortex core, so 
that it is better to dismiss the center of the circular structures totally, 
concentrating instead on flat nano-rings, i.e. with an inner radius Ri 
and an outer radius i?2- In such nano-rings, the stability of the vortex 
state is enhanced, and concerning the switching of the vorticity, they 
have similar properties as circular ones, i.e. with Ri = 0. 

For in-plane flelds, the fleld strength need only to be > 200 Oer- 
steds, but the time scales of vorticity-switching are > 1 ns; the vortex 
center is typically unstable, i.e. the topology of the state is changed 
from a 'deformed vortex state' to a 'deformed in-plane single-domain 
state', if in-plane flelds of the 'critical duration' arc applied. So the 
case of 'in-plane vortex switching flelds' should be dismissed rightaway 
for circular dots. 
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Figures and Figure Captions 



(^(r) □ OOMMF 




X [ nm] 



Fig. la: The 'profile function' az{r) as given by our analytical ansatz 
is compared with results of the direct numerical simulation by the 
OOMMF code. Material parameters: Permalloy. The radius of the 
circular dot considered is 175.5 nm, the height is 6 nm. 
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Fig. lb: The effective radius -Refr of the vortex core is presented over 
the reduced thickness := p- of the circular dot, in units of the 
magnetic exchange length 1^- Material parameters: Py, as in Fig. la. 
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Figure 2: The critical thickness tc{R), i.e. the hne which separates 
the stabihty ranges of the vortex state from that of the in-plane single- 
domain state, is presented as a function of the radius R of the circu- 
lar dot. Parameters of Permalloy material are used. The solid line 
describes our results, as given in the text; the dashed line presents 
the simplified approximation tcfl(R), also described in the text, and 
the circles describes the results of a numerical calculation with the 
OOMMF code. (Besides, the additional phase transition line separat- 
ing the out-of-plane single-domain state from the in-plane single do- 
main state would in the present plot almost coincide with the abscissa 
axis.) 
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Figure 3: The total energy density / ( the sohd hne), and the par- 
tial energy densities of the dipolar energy (squares), exchange energy 
(diamonds) and Zeeman energy (triangles)) for the vortex state in a 
circular permalloy dot of radius R = 150 nm (note the vertical line) 
and thickness th = 4: nm are presented as a function of the shift A of 
the vortex core; for more details see the text 
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Figure 4: For different values of the external in-plane Zeeman field the 
density /tot of the total energy of Fig. 3 is presented against the shift 
A of the vortex core. One can see that for strong enough Zeeman fields 
{Hy^^ > 190 Oe) there is no longer a local minimum, i.e. the vortex 
vanishes (A — »• oo), which means a transition to a single-domain state. 
The vertical line denotes the radius R of the circular dot. 
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Figure 5: For Py dots of thickness th = 15 nm and 8.3 nm, respec- 
tively, experimental results, and our theoretical results, for the in-plane 
annihilation fields of the vortex state are presented over the double re- 
ciprocal aspect ratio 2R/th of the dot. 
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Figure 6: The critical time tc for dynamic annihilation of the vortex 
state by in-plane fields Hy^^ is presented against H^^^ for different values 
of the Gilbert damping ac- For further details see the text. 
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Figure 7: The line separating the dynamic stability and instability 
regions for dynamic annihilation of the vortex core by in-plane fields, 
as in the preceding figure, is presented in an (ac Hy^^) plane, where 
ac is Gilbert's damping parameter and H^^^ the in-plane 'switching 
field'. 
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Figure 8: This 'movie sequence' presents the annihilation of a vortex 
state by an in-plane field. This is a rather slow process, which needs 
at least 1.5 ns, although the in- plane field is as strong as 300 Oe . 
Gilbert's damping parameter has been assumed to be aa — 0.02 



26 



100 




[kOe] 

Figure 9: This plot presents the typical 'switching time' ti/2 of a 
vortex state by a pulsed perpendicular field Hf^. This is a very fast 
process, which needs typically only ~ 40 ps, but the necessary strength 
of the field is very large, ~ 5000 Oe. Gilbert's damping parameter has 
been assumed to be ac = 0.02. "Theory" means the simple result 
ti/2 — '^/{\lo\Hf^); the other results are from OOMMF-simulations. 
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Figure 10: This plot sliows details of the fast switching of the vor- 
ticity of our flat circular Permalloy dot of radius R = 150 nm and 
thickness th = 4: nm by a pulsed vertical field of Hf^^ = 5000 Oe. The 
energy densities / of the dipolar energy ('stray-field') and of the ex- 
change energy are presented over the time for the first 50 ps after the 
start of the perpendicular field pulse, which lasts from t = to t = 42 
ps. 
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Figure 11: This plot sliows tiie 'wiggling' of the small out-of-plane 
component of the (spatially averaged) magnetization of the circular 
Permalloy dot; comes to rest only after several ns, although the 
switching of the vorticity needs essentially only 42 ps, see the preceding 
figure. Again we have assumed = 0.02. 



29 



Ops 

































*- 


\ 






i 




>^ 


\ 




1 

* 


i 
t 


T T T 


^\ \ 

■ * M 


\ 
• 


\ 
• 


\ 


\ 
\ 


i ^ \' 
\ 


>^ t f 
*^ ' f 


> 


f 




\ 



























6( 


5] 


ps 

..*s<. ... 


\ 




JT- 

ji' *- ^- 




\ 


^ t 


\ \ 
\ 1 










\ t 


f ^ 


i; 

t: 


If ' 


i 4 






f # 


1 'f 






\ \ 






• f 


r f 




V 


\ \ 












V 










• t 






•4 


-It 







170 ps 





4- -4- Vv. 






/ 
; 


/ . M- \ ^ 
tit 






* t f f t 


T 


V V \ ^ 




^ \ 



220 ps 500 ps 800 ps 









4-- 








































*- 


4 








* 
























* 










r 


«- 






^ 










^ 








\ 














If- 


If 










:t: 




s 










if 


-\ \ \ 








; / / 




\ 




* 


*' 




*- 




f 


-!•- 

-.1.- 


* 






t 
r 




i 


I *' 


-^ t t 


1 


t 


\ 








/ 






i 






> 


* 






f 

1 


I 


t 


I 


^ \ 


■^f f f 






\ 










t 




i 




\ 


\ 










f 


1 


\ 




\ 






/ 


1 
* 


\ 




■^-*/< f * 


4 

1 






\ 


iV 


'v * 














\ 












\ 






; 






\ 


^( 


^* 


-* 








• 






\ 








* 










/ 






hi 






b 














* 




--» --»■ ^ .-r 






































































1" 


















-•-ir... 

















Figure 12: Here white (black) background colours mean positive (neg- 
ative) values of a^; so this 'movie sequence' shows that under the 
conditions described in the text the vorticity is not switched at all, 
whereas the central polarity changes from +1 to -1 at t ~ 60 ps, re- 
mains there until i ~ 170 ps, but returns to +1 for t > 220 ps. The 
applied field pulses were H^^^ = 1 kOe between and 60 ps, and si- 
multaneusly H^^^ = 150 Oe between and 30 ps. As before, the radius 
and the thickness of our Py dot were 150 and 4 nm, respectively, and 
the Gilbert damping was ao — 0.02. Note that for 60 ps and 220 ps 
one is no longer dealing with a vortex state, but rather with a spiral. 
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Figure 13: For flat 'nano-rings' made of Py material with fixed outer 
radius R2 = 150 nm and thickness = 4 nm, the total energy densities 
/tot of the vortex state and the homogeneous in-plane state ('single- 
domain state') are presented against the inner radius Ri. 
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Figure 14: The switching time ti/2 of the vorticity of a fiat circular 
Permalloy dot (thickness th = 4 nm, radius R = 150 nm) and of a 
corresponding 'nano-ring' with outer radius R2 = R and inner radius 
Ri = R/3 is plotted against the vertical switching field Hf^^. 
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Figure Captions 



Fig. la: The 'profile function' a^ir) as given by our analytical ansatz 
is compared with results of the direct numerical simulation by the 
OOMMF code. Material parameters: Permalloy. The radius of the 
circular dot considered is 175.5 nm, the height is 6 nm. 
Fig. lb: The effective radius Res of the vortex core is presented over 
the reduced thickness Th := -p- of the circular dot, in units of the 

I'm 

magnetic exchange length Material parameters: Py, as in Fig. la. 
Figure 2: The critical thickness tc{R), i.e. the line which separates 
the stability ranges of the vortex state from that of the in-plane single- 
domain state, is presented as a function of the radius R of the circu- 
lar dot. Parameters of Permalloy material are used. The solid line 
describes our results, as given in the text; the dashed line presents 
the simplified approximation tcfi{R), also described in the text, and 
the circles describes the results of a numerical calculation with the 
OOMMF code. (Besides, the additional phase transition line separat- 
ing the out-of-plane single-domain state from the in-plane single do- 
main state would in the present plot almost coincide with the abscissa 
axis.) 

Figure 3: The total energy density / ( the solid line), and the par- 
tial energy densities of the dipolar energy (squares), exchange energy 
(diamonds) and Zeeman energy (triangles)) for the vortex state in a 
circular permalloy dot of radius R — 150 nm (note the vertical line) 
and thickness th — 4: nm are presented as a function of the shift A of 
the vortex core; for more details see the text 

Figure 4: For different values of the external in-plane Zeeman field the 
density /tot of the total energy of Fig. 3 is presented against the shift 
A of the vortex core. One can see that for strong enough Zeeman fields 
(Hy^ > 190 Oe) there is no longer a local minimum, i.e. the vortex 
vanishes (A oo), which means a transition to a single-domain state. 
The vertical line denotes the radius R of the circular dot. 
Figure 5: For Py dots of thickness — 15 nm and 8.3 nm, respec- 
tively, experimental results, and our theoretical results, for the in-plane 
annihilation fields of the vortex state are presented over the double re- 
ciprocal aspect ratio 2R/th of the dot. 
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Figure 6: The critical time tc for dynamic annihilation of the vortex 
state by in-plane fields Hy^*" is presented against for different values 
of the Gilbert damping ac- For further details see the text. 
Figure 7: The line separating the dynamic stability and instability 
regions for dynamic annihilation of the vortex core by in-plane fields, 
as in the preceding figure, is presented in an (ckg, H"^^) plane, where 
aa is Gilbert's damping parameter and Hy^^ the in-plane 'switching 
field'. 

Figure 8: This 'movie sequence' presents the annihilation of a vortex 
state by an in-plane field. This is a rather slow process, which needs 
at least 1.5 ns, although the in- plane field is as strong as 300 Oe . 
Gilbert's damping parameter has been assumed to be aa — 0.02 
Figure 9: This plot presents the typical 'switching time' ti/2 of a 
vortex state by a pulsed perpendicular field H'^'^. This is a very fast 
process, which needs typically only ~ 40 ps, but the necessary strength 
of the field is very large, ~ 5000 Oe. Gilbert's damping parameter has 
been assumed to be = 0.02. "Theory" means the simple result 
ti/2 = 7r/(|7o|//^^*); the other results are from OOMMF-simulations. 
Figure 10: This plot shows details of the fast switching of the vor- 
ticity of our fiat circular Permalloy dot of radius i? 150 nm and 
thickness t/, = 4 nm by a pulsed vertical field of Hf^^ = 5000 Oe. The 
energy densities / of the dipolar energy ('stray-field') and of the ex- 
change energy are presented over the time for the first 50 ps after the 
start of the perpendicular field pulse, which lasts from t — to t — 4:2 
ps. 

Figure 11: This plot shows the 'wiggling' of the small out-of-plane 

component of the (spatially averaged) magnetization of the circular 
Permalloy dot; comes to rest only after several ns, although the 
switching of the vorticity needs essentially only 42 ps, see the preceding 
figure. Again we have assumed ac; = 0.02. 

Figure 12: Here white (black) background colours mean positive (neg- 
ative) values of a^; so this 'movie sequence' shows that under the 
conditions described in the text the vorticity is not switched at all, 
whereas the central polarity changes from -|-1 to -1 at t 60 ps, re- 
mains there until t ^ 170 ps, but returns to +1 for t > 220 ps. The 
applied field pulses were Hf^^ = 1 kOe between and 60 ps, and si- 
multaneusly Hy^^ = 150 Oe between and 30 ps. As before, the radius 
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and the thickness of our Py dot were 150 and 4 nm, respectively, and 
the Gilbert damping was ac = 0.02. Note that for 60 ps and 220 ps 
one is no longer dealing with a vortex state, but rather with a spiral. 
Figure 13: For flat 'nano-rings' made of Py material with fixed outer 
radius i?2 = 150 nm and thickness = 4 nm, the total energy densities 
/tot of the vortex state and the homogeneous in-plane state ('single- 
domain state') are presented against the inner radius Ri. 
Figure 14: The switching time ti/2 of the vorticity of a fiat circular 
Permalloy dot (thickness th = 4 nm, radius R = 150 nm) and of a 
corresponding 'nano-ring' with outer radius R2 = R and inner radius 
Ri — R/S is plotted against the vertical switching field Hf^^. 
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